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1.  SC  vehicle  as  a result  of  tendencies  of  high  speed  vehicle  development 


In  practice,  the  vehicle  development  is  stipulated  by  a useful  effect  that  is  attained  upon  the  vehicle 
operation. 

The  useful  effect  for  transport  ships  consists  in  cargo  transportation.  Thus,  the  ship  efficiency  as 
vehicle  efficiency  will  be  determined  by  commercial  effectiveness  [22] 
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where  P is  the  payload,  V is  the  velocity,  N is  the  power. 

The  velocity  and  range  are  determining  parameters  for  objects  used  for  military  and  rescue  purposes. 
In  this  case  one  can  see  that  the  hydrodynamic  drag  reduction  positively  influences  on  all  main 
characteristics  and  increases  the  commercial  effectiveness  coefficient  owing  to  consumed  power 
decrease  or  increase  of  the  velocity  and  range. 

For  surface  vehicles,  we  clearly  see  a tendency  of  decreasing  the  waterline  area  and  wetted  surface  of 
the  hull  due  to  dynamic  means  of  maintenance  of  the  hull  on  the  water  surface.  There  are  ships  with 
small  waterline  area  SWATHS,  hydrofoil  ships,  hovercrafts,  ekranoplanes  and  hydroplanes. 

In  this  case  the  energy  consumption  to  realize  the  dynamic  means  of  maintaining  the  hull  on  the  water 
surface  may  be  related  to  the  controlling  effect  for  the  drag  reduction. 

One  can  see  in  Fig.  1 that  the  velocity  and  specific  power  of  the  vehicles  increases  at  exiting  the  hull 
from  the  water. 

The  same  way  of  reduction  of  the  wetted  surface  of  the  hull  due  to  developed  artificial  or  vapor 
cavitation  - supercavitation  [10,  11,  15]  may  be  used  for  underwater  motion. 

Calculations  shows  it  is  possible  to  obtain  a twenty-fold  advantage  in  the  hydrodynamic  drag  when 
the  velocity  is  equal  to  100  m/sec,  and  cavitation  number  o = 0.01,  and  1000-fold  advantage  may  be 
obtained  when  the  cavitation  number  a = lO^1  [11]. 


Po  ~ P 

v2  r. 


where  /70  is  the  hydrostatic  pressure,  pc  is  the  cavity  pressure,  which  is  equal  to  the  saturated  water 
vapor  pressure  pc  = pv  at  vapor  cavitation. 


Paper  presented  at  the  RTOAVT  Lecture  Series  on  “Supercavitating  Flows  ",  held  at  the  von  Kdrmdn 
Institute  (VKI)  in  Brussels,  Belgium,  12-16  February  2001,  and  published  in  RTO  EN-010, 
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Comparing  the  efficiency  of  the  different  methods  of  the  hydrodynamic  drag  reduction  with  efficiency 
of  the  artificial  and  vapor  supercavitation,  one  can  see  its  advantage  and,  hence,  perspective  of  its  use 
[22],  Moreover,  a vapor  layer  between  the  water  and  hull  unites  many  methods  of  the  drag  reduction: 

• Reduction  of  water  adhesion  to  the  surface  or  slippage; 

• Effect  of  the  movable  surface  (the  gas  flow  velocity  in  the  clearance  may  be  equal  to  the  water  flow 
velocity); 

• Change  of  the  physical  constants  of  BL  (here,  the  medium  density  reduction); 

• Gas  blowing  into  the  boundary  layer  (BL). 

Here,  - in  the  case  of  the  artificial  cavitation,  the  gas  is  blown  through  the  solid  surface  of  the  object. 

- in  the  case  of  vapor  supercavitation,  the  gas  arrives  from  the  supercavity  surface. 

However,  use  of  the  supercavitation  flow  schemes  requires  respectively  high  specific  power  of  the 
vehicle.  So,  when  the  velocity  is  equal  to  100  m/sec,  and  cavitation  numbers  a = 0.05  h-  0.01,  specific 

available  power  will  be  of  order  -^-  = 2000 4000  H.p./t  This  is  compared  with  specific  power  of  the 

planing  ships  in  the  universal  diagram  of  Gabrielly-Karman  [22]  (Fig.  2) . 


2.  Features  of  supercavitation 

First  of  all,  the  supercavitation  flow  regime  gives  limitations  for  shape  and  dimensions  of  the  hull  of 
an  underwater  vehicle.  They  must  correspond  to  the  calculated  supercavity  shape.  Therefore, 
formalization  of  the  hull  geometry  will  be  performed  by  comparing  the  vehicle  volume  and  surface  with 
the  supercavity  volume  and  surface  or  with  shapes  of  hulls  of  Reichhardt  series  having  the  constant 
pressure  [19], 

In  this  case  dimensions  of  the  vehicle  having  a shape  of  a body  of  revolution  will  be  restricted  by 
limiting  dimensions  of  the  supercavity.  They  are  determined  by  known  formulae  for  main  supercavity 
dimensions  [1] 


where  cavitation  number  cr  = 2{p0  — pc)/  pV 2 is  a main  parameter  of  the  supercavitation  flows, 
p(),  p.  are  pressures  in  the  stream  and  in  the  cavity,  respectively,  and  cx  = cx0  (l  + O')  is  the  coefficient 
of  this  cavitator  in  the  form  of  a disk  with  diameter  Dn , cr0  is  the  drag  coefficient  when  <7  = 0. 

We  should  note  that  the  supercavitation  motion  regime  requires  enough  high  motion  velocity  that  is 
sufficient  to  form  a supercavity  where  the  vehicle  hull  is  inscribed  completely.  This  is  explained  by 
structure  of  the  formulae  for  main  supercavity  dimensions  (1): 
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The  object  calculated  for  given  regime  and  inscribed  into  the  calculated  supercavity  will  need  the 
attainment  of  this  calculated  velocity  in  the  unseparated  flow  regime.  During  acceleration  some  difficult 
problems  arise.  These  are: 

- energetic  problem  connected  with  overcoming 

the  high  hydrodynamic  drag  of  the  vehicle  in  unseparated  flow; 

- problem  on  ensuring  the  motion  stability  in 
non-calculated  flow  regime. 

The  second  essential  feature  of  the  supercavitation  motion  is  a necessity  to  take  into  account  the 
energy  expenditure,  volume  and  mass  to  ensure  artificial  regime  of  supercavitation,  i.e.  gas  blowing  into 
the  cavity.  This  is  inevitable  upon  the  object  acceleration  up  to  the  calculated  velocity  and  maintaining 
the  supercavity  in  sub-vapor  supercavitation  regimes  when  the  pressure  in  the  cavity  is  higher  that  the 
pressure  of  saturated  water  vapors.  The  necessity  of  the  gas-supply  into  the  cavity  disappears  at 
attainment  of  vapor  supercavitation  regimes  when  the  supercavity  is  filled  by  water  vapor. 

The  problem  of  interaction  between  the  propeller  and  supercavity  is  sensitive  to  the  velocity  and 
pressure  perturbations  caused  by  working  propeller.  This  may  result  in  additional  washing  the  hull  and 
increasing  the  gas-supply  rate. 

In  the  general  case  the  supercavitating  vehicle  motion  may  be  considered  as  one  consisting  of  three 
phases  : 

- phase  1 - the  vehicle  acceleration  up  to  the  calculated  velocity; 

- phase  2 - motion  with  constant  velocity; 

- phase  3 - motion  on  inertia. 

The  known  applications  of  the  supercavitation  motion  are  used  to  attain  the  maximal  motion  range  in 
phase  2 or  in  phase  3.  In  this  case  phase  1 is  not  essential  to  pass  the  distance,  but  has  decisive 
significance  in  the  propulsive  sense,  since  requires  development  of  the  highest  power  exceeding  multiply 
the  power  necessary  for  motion  with  constant  calculated  velocity. 


3.  Energy  of  SC  vehicles  in  stationary  motion 


Energy  of  SC  vehicles  in  the  stationary  motion  will  be  determined  by  their  hydrodynamic  drag  in 
steady  flow.  Since  the  vehicle  hull  will  place  in  a vapor  or  gas  cavity,  the  main  component  of  the 
hydrodynamic  drag  is  the  cavitator  drag.  It  is  determined  by  formula 
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pv2 

2 


r 


(3) 


where  p is  the  fluid  density,  V is  the  motion  velocity,  cx  is  the  cavitator  drag  coefficient  (for  disk 

kD] 

cx  = 0.82(1+  (7) , where  <3  is  the  cavitation  number),  Sn  = — is  the  cavitator  area. 

Additional  hydrodynamic  drags  can  arise  owing  to  using  the  schemes  of  supercavitation  flow  with 
partial  or  periodic  washing  the  tail  part  of  the  object.  Fig.  3 shows  possible  schemes  of  supercavitation 
flow  around  SCO: 

a)  scheme  of  flow  with  forming  two  cavities; 

b)  stationary  planing  along  the  inner  cavity  surface; 
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c)  periodic  impact  interaction  with  the  cavity  walls; 

d)  aerodynamic  interaction  with  vapor-splash  medium  of  the  cavity. 

In  the  general  case  the  total  hydrodynamic  drag  of  SCO  may  be  represented  with  taking  account  of 
the  cavitator  drag  ( 3)  in  the  form 


K=Rx„+RxF,  (4) 

where  RxF  is  the  additional  hydrodynamic  drag  of  friction  upon  washing  the  hull.  It  may  be  equal 
Rxn  = 0.0 \Rxn  at  the  flow  scheme  d)  (Fig.  3)  and  Rxn  = 0.5 Rxn  at  the  flow  scheme  a)  (Fig.  3). 

When  V = Const , the  propeller  thrust  T must  equalize  the  drag  (4) 


T = RX- 

If  the  propeller  is  one  of  rocket  type,  its  thrust  will  be  equal 

T = mfQf, 


where  m f = dm  / dt  is  the  mass  rate  of  the  fuel.  O . is  the  specific  impulse  of  the  fuel. 
For  steady  motion  when  mf  = Const, 


T 


Qf 


(5) 


(6) 


(7) 


Using  the  formulae  (3)  - (7),  we  obtain 


V = 


(8) 


Introducing  the  passed  distance  S in  the  form 


S = Vt. 


(9) 


and  define  the  volume  due  to  the  fuel 


Wf  = mf  / p f = const , 


(10) 
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we  obtain  for  the  passed  distance  that 


5 = 


~~  C^PfQft , 


(11) 


2 Wf 

where  — = C,  = Const , p . is  the  fuel  density. 

\cxpSn 


One  can  see  in  the  expression  (11)  that  to  attain  the  maximal  range  S it  is  necessary  to  apply  fuels 
with  maximal  specific  impulse  Q f and  density  p f . 

The  vehicle  velocity  is  obtained  from  expressions  (3)  and  (5) 


V = 


\ 


(12) 


We  can  see  from  here  that  the  velocity  is  proportional  to  ratio  (T  / Sn  )V~ . 

We  can  obtain  the  estimation  of  maximal  range  from  the  assumption  that  whole  cavity  volume  is  filled 
by  the  vehicle  hull,  and  the  hull  is  filled  by  fuel  with  density  p f and  specific  impulse  Q f . 

As  since  the  cavities  have  a shape  close  to  an  ellipsoid  with  known  diameter  and  length  (1),  it  is 
possible  to  calculate  the  supercavitry  volume  Wc  with  taking  account  of  the  relation  Lc  = XDC  (1) 


we  =0.71tp  R;  lc. 


(13) 


The  SCO  drag  will  be  equal  [1] 


R 


X 


kScAP , 


(14) 


where  = nR;  ,s  the  cavrty  m.d-secfion  area,  AP  = Po  - = CT—  is  the  pressure  difference, 

k = 0.96 . 


The  fuel  mass  in  the  cavity  volume 


mf=pfWc=0Jl7iR2cLcpf. 


(15) 
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Time,  when  the  propeller  is  maintaining  the  necessary  thrust,  or  the  motion  time  is 


(16) 


where  Oy  = m , Qf  is  the  total  impulse  of  the  fuel. 


Hence,  the  maximal  range  of  the  motion  S will  be  equal 


S = 


tV  = \ A— 

p V o 


(17) 


or  in  the  dimensionless  form 


S=  — = (18) 

4 p v 

If  we  assume  that  metal  fuel  of  high  efficiency  is  used  ( p f / p = 2.7;  Qf  IV  = 30;  <J  = 0.01;  Lc  = 10  m), 
then  the  maximal  range  will  be  equal  to 


S = 1 .4  • 1 00  • 30  • 2.7  = 1 1 340 

or  in  dimensional  form  S = 1 13400  m. 

When  the  solid  rocket  fuel  on  the  basis  of  nitrocellulose  and  nitroglycerin  - JPN  with  density 
p f = 1 .6;  Qf  / V = 1 1 is  used,  then  the  maximal  range  will  be  equal 


S = 1.4-1 00  -1.6-11  = 2464 , 

S = SLC  = 24640  m. 

Calculation  of  maximally  attainable  range  is  represented  in  graphic  form  in  Fig.  4.  Here,  the  range  S is 
represented  as  a dependence  on  velocity  V and  aspect  ratio  of  the  hull  (cavity)  X = LC  / Dc  for  specific 
impulses  of  the  fuels  Qf  = 230  sec  and  Qf  = 600  sec,  depth  H = 10  m,  Lc  = 10  m.  If  it  is  assumed 
that  the  strengthening  (design)  restriction  of  the  hull  aspect  ratio  is  X < 20,  then  the  maximal  restriction 
on  the  motion  range  of  SCO  is  Sm,dX  = 70  km  and  on  the  velocity  Fmax  = 200  m/sec. 


3.1.  Underwater  vehicle  propellers 

Underwater  vehicle  propulsion  is  intended  to  transform  chemical  energy  of  the  fuel  into  thermal 
energy.  The  latter  transforms  into  mechanical  energy  - work  for  overcoming  the  hydrodynamic  drag  Rx 
(3). 
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When  the  motion  is  uniform,  the  thrust  T of  the  propeller  must  equalize  the  hydrodynamic  drag  Rx 

T = RX.  (19) 

In  this  case  the  useful  power  of  the  underwater  vehicle  propulsion  will  be  equal  to 

N = TV.  (20) 

The  propulsion  effectiveness  is  estimated  by  its  total  efficiency  Tj  representing  the  ratio  of  useful  power 
N and  all  consumed  power 

TV 

*7  = — ’ (21) 

I q 

where  I is  the  mechanical  equivalent  of  heat;  q is  the  heat  rate  per  time  unit. 

The  coefficient  r)  characterizes  the  propulsion  perfection  and  is  dependent  of  the  fuel  characteristics. 
In  practice,  the  specific  rate  of  the  fuel  is  of  importance.  It  is  quantity  of  the  fuel  consumed  per  time  unit 
at  unit  power.  It  is  usually  measured  in  g/kilo-Watthour.  Since  q = qmf,  where  q is  the  caloric 

ability  of  the  fuel;  m , is  the  fuel  rate  per  a second. 

Using  (21),  we  obtain  the  relation  for  the  specific  rate  of  the  fuel  and  efficiency  Tj 


mf  _ l 

~TV~Tqrj' 


(22) 


One  can  see  from  (22)  that  the  specific  rate  of  the  fuel  q , depends  both  on  the  propulsion  characteristics 
Tj  and  on  the  fuel  characteristics  q . For  diesels 

qf  =215  g/kW  • h,  for  turbines  qf=  340  g/kW  • h. 

We  well  know  specific  powers  of  different  propulsions  used  in  the  underwater  vehicles  [4].  There  are 

diesel  propulsions  - 45 + 90  W/kg  |10  30  W/  m3 1 

electro-accumulators  -60+100  W/kg  | 200  | 

atomic  propulsions  - 25  + 200  W/kg  | 11.3  | 

high-speed  diesels  - 1.33  kW/kg  | 224  kW/ nr  j 

gas  turbines  - 5.7  kW/kg  | 1700kW/m3| 

It  is  known  as  well  that  the  turbine  or  piston  propulsion  volume  does  not  exceed  50%  of  the  vehicle 
volume,  whereas  the  fuel  part  in  rockets  exceeds  90%. 
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It  is  essential  that  the  necessary  specific  power  for  supercavitating  vehicles  is  calculated  by  formula 


RXV  RXV  nnxopV2 

pv  = — — = = 0.71— 

Wc  0.61k  R2LC  Lc 


(23) 


When  o = 0.01,  Lc  = 1 0 m,  pv  = 715  kW/m3  at  V = 100  m/sec;  pv=5720  kW/m3  atF=200 
m/sec. 

These  values  exceed  on  2 - 3 orders  the  analogous  values  for  diesel,  electric,  atomic  propulsions. 
Therefore,  the  gas-turbine  and  jet  propulsions  with  the  metal  fuel  and  using  outboard  water  as  a oxidizer 
and  a cooler  of  the  combustion  products  are  considered  perspective  ones  for  the  supercavitating  vehicles 
at  increased  motion  velocity.  Such  propulsions  were  investigated  experimentally,  and  their  high 
efficiency  was  confirmed  [19], 

At  that,  the  increase  of  the  velocity  of  SCO  motion  compared  with  vehicles  at  unseparated  flow  up  to 
100  - 200  m/sec  do  attractive  the  rocket  schemes  with  solid  fuel  in  view  of  the  propulsion 
efficiency  increase 

(Fig.  5) 


y,-y 

vj 


(24) 


Fig.  5 represents  a graph  of  dependence  of  the  propulsive  efficiency  for  the  jet  propeller  with  solid 
fuel  on  the  motion  velocity  at  Vj  = 2040  m/sec. 

Reaction  ,4/ + water  gives  high  temperature  10800K.  Therefore,  water  must  be  supplied  into  the 
propeller  in  plenty.  In  this  case  the  propulsion  temperature  reduction  will  occur  according  to  the 
theoretical  dependence  on  water  gram-molecules  n (Fig.  6a). 

It  is  well  known  that  to  accelerate  the  reaction  it  is  necessary  to  fluidize  the  metal  and  to  transform 
the  water  into  vapor.  The  calculation  of  the  thermal  processes  is  performed  at  condition  that  the 
temperature  of  the  reaction  products  was  reduced  to  1250K.  A scheme  of  the  propulsion  working  on  fuel 
A l with  using  the  outboard  water  is  shown  in  Fig.  6b.  In  the  presented  propulsion,  the  heat  from 
combustion  chamber  is  transmitted  to  aluminum  to  form  a melt  at  temperature  about  950K  and  to  water 
to  form  vapor.  The  combustion  products  arrive  on  the  turbine  and  after  the  turbine  to  outboard.  In  this 
case  the  propulsion  efficiency  will  depend  essentially  on  the  excess  pressure  in  the  combustion  chamber 
on  the  different  depths.  This  is  shown  in  graph  in  Fig.  7. 

The  graph  in  Fig.  7 gives  the  dependence  of  the  specific  output  power  of  the  propulsion  with  the  fuel 
- aluminum  + sea- water  for  different  depths  H=  305  + 1525  m.  For  supercavitation  motion  regimes,  the 
vapor  exhaust  will  occur  into  the  supercavity  at  pressure  P = 0.02  kg/cm3.  This  corresponds  to  the 
saturated  vapor  pressure  at  temperature  1 5°C.  This  pressure  corresponds  to  the  air  pressure  on  height  23 
km  above  sea  level  and  ensures  ideal  conditions  for  work  of  propulsions  with  opened  cycle  independently 
on  the  object  motion  depth. 

A scheme  of  the  turbo-rotor  propeller  for  supercavitating  vehicle  was  considered  to  increase  its 
efficiency  at  the  Institute  for  Fluid  Mechanics.  Propellers  with  one  or  two  turbo-rotors  (Fig.  8,  a,b,c)  were 
tested  there.  It  is  shown  in  the  scheme  that  the  gas  turbine  and  screw  with  three  blades  working  in 
supercavitation  regime  were  combined  on  one  disk  of  the  rotor.  There  also  are  shown  section  of  the 
turbine  blades  and  supercavitating  screw  sections. 
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Photographs  in  Fig.  9 a,b,c  show  models  of  the  turbo-rotors  in  the  hydro-tunnel  and  a picture  of  flow 
around  the  working  turbo-rotors  in  the  hydro-tunnel  when  the  free  stream  velocity  V=  10  m/sec  at  using 
the  gas  turbine  with  compressed  air  and  at  pressure  to  6 kg/cm2. 

Preliminary  experiments  with  turbo-rotor  propellers  showed  real  possibility  of  the  propeller  thrust 
increase  on  20  % compared  to  the  rocket  scheme,  although  theoretical  estimations  point  to  possibility 
of  double  increase  of  the  thrust.  Therefore,  possibility  of  attainment  of  optimal  conditions  for  work  of  the 
turbine  and  supercavitating  screw  still  will  have  to  be  realized. 


3.2  Fuel  for  SCO 


We  can  see  from  the  previous  analysis  that  the  fuel  for  high-speed  SCO  must  satisfy  the  heightened 
requirements  to  the  specific  characteristics  Q f and  p , . The  specific  characteristic  - power/volume  is 
the  more  important  than  the  power/mass. 

Work  [21]  considers  schemes  of  turbo-prop  of  the  hydro-jet  propeller  working  on  unitary  solid  fuel 
(Fig.  10,  a),  on  two-component  fuel  (Fig.  10,  b)  and  on  hydro-reacting  fuel  on  the  basis  of  reaction 
between  the  magnesium  and  sea-water  (Fig.  10,  c). 


Calculations  showed  that  ratios  of  specific  rates  of  the  fuel  transported  on  the  board 
wide  velocity  range  20  4-  100  m/sec  have  constant  values 


m f 

rhfa=—  in 


m ( schem  b)  m c,,  (schem  c) 

— = 0.42;  -i - = 0.37. 

m fa  ( schem  a)  m )a  ( schem  a) 

The  relations  show  that  the  hydro-reacting  fuel  on  the  vehicle  board  must  be  in  1.15  times  less  than 
the  two-component  one  (scheme  b)  and  in  2.7  times  less  than  the  solid  unitary  fuel  (scheme  a). 

In  rocket  systems,  where  the  fuel  mass  predominates,  and  weight  and  volume  of  the  nozzle  and 
combustion  chamber  are  comparatively  small,  the  propulsive  characteristics  of  the  system  will  be 
determined  by  the  fuel  characteristics  - its  specific  impulse  Q/  and  density  p , . 

The  investigations  showed  that  the  chemical  fuels  using  metals  and  sea-water  as  oxidizer,  ballast  and 
cooler  are  the  most  effective  for  underwater  vehicles.  This  is  stipulated  that  the  metals  ensures  maximal 
energy  at  reaction  with  water  [ 19,  21] 


Al  + sea-water  — > 3 1 43  kW  hour/m3  (1.18  kW  hour/kg), 
Mg  + sea-water  — > 1891  kW  hour/m3  (1.05  kW  hour/kg), 
Li  + sea-water  — > 1092  kW  hour/m  3 (2.09  kW  hour/kg). 


Reactions  and  main  parameters  are  presented  in  Table  [ 19,21]  Fig.  11. 

In  practice,  Al  is  the  best  fuel,  since  it  is  proof  and  cheap.  Aluminum  melting  temperature  is  950K, 
and  the  combustion  temperature  in  reaction  with  H20  reaches  10800K. 

Perspective  of  the  fuels  for  the  underwater  vehicles  is  frequently  estimated  by  theoretical  value  of 
specific  work  of  the  jet  flow  (kW-h/m  J ).  It  is  determined  at  fixed  value  of  the  balanced  expansion  and 
quantity  of  the  injected  sea-water  sufficient  to  decrease  the  combustion  product  temperature  to  1250K. 
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Dependence  of  the  approximate  theoretical  values  of  the  specific  power  with  respect  to  mass  p,  and 
volume  p2  for  Al  at  reaction  with  water 


2AI  + n ■ H20,  (25) 

where  n is  the  number  of  gram-molecules  upon  expansion  of  the  combustion  products  from  206.9  to 
10.14  kN/m2 , is  presented  in  Fig.  12. 


4.  Calculation  of  supercavitating  model  motion  range  at  inertial  motion 


Drag  of  a supercavitating  model  completely  enveloped  by  a cavity  is  equal  to  a cavitator  drag  and 
calculated  by  formula 


R = 


PV2 


r 9 

^xn^n  9 


(26) 


where  p is  the  fluid  density;  V is  the  model  velocity;  cxn  is  the  cavitator  drag  coefficient  (for  disk  cxn  = 
0.82  ( 1 + a ),  where  o is  the  cavitation  number);  Sn  is  the  cavitator  area. 


When  cavitation  number  is  small,  it  is  possible  to  consider  cxa  as  constant  one.  Then  the  equation  of 
motion  of  the  model  with  mass  m has  the  form 


dV 

dt 


-kV 


2 


(27) 


where  k ■ 


P cvA 
2m 


Its  solution  gives  dependence  of  motion  time  on  velocity  in  the  form 


t(V)  = - 
k 


V~K 


and  dependence  of  passed  range  S on  the  velocity  in  the  form 


1 V 

S(V)  = - In-*-  , 
k V 


(28) 


(29) 


where  V0  is  the  starting  velocity  of  the  model. 
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Hence,  the  model  moves  infinitely  long  and  passes  infinitely  great  distance.  It  is  necessary  to  define 
criterion  determining  the  finite  point  of  the  path  to  have  sense  of  the  range  calculation. 

Such  criteria  can  be: 

- the  velocity  reduction  to  certain  value; 

- decrease  of  the  model  momentum  to  certain  value; 

- disappearance  of  clearance  between  the  model  and  cavity  walls. 

In  essence,  the  second  and  third  criteria  are  distinctive  from  the  first  one  only  by  way  of  the  final 
velocity  definition.  Thus,  it  is  possible  to  consider  the  final  velocity  Ve  definition  as  the  most  universal 

way  of  determination  of  the  motion  range  Se . 

The  range  depends  on  the  model  mass  m and  cavitator  radius  rp  determining  the  drag  for  given 
velocities  V0  and  Ve. 

As  the  model  must  be  located  in  the  cavity,  its  maximal  possible  mass  depends  on  the  material  density 
pb  , model  length  L and  cross  cavity  dimensions  determined  by  the  cavitator  radius. 

Thus,  the  range  is  determined  by 

- material  density; 

- model  length; 

- cavitator  radius. 


Calculations  necessary  to  determine  the  range  can  be  executed  due  to  the  program  SCAV  [23].  The 
model  shape  is  approximated  by  three  truncated  cones  (Fig.  13,  a)  in  this  program.  The  program 
computes: 

- the  model  characteristics:  volume,  mass  center  location  and  moment  of  inertia; 

- the  shape  and  dimensions  of  stationary  supercavity  past  a disk  or  a conic  cavitator  at  given  velocity 
and  depth  of  the  motion; 

- the  clearances  between  the  model  and  cavity  boundary; 

- the  value  of  impact  load  acting  to  the  model  at  water  entry; 

- the  speed  reduction  of  the  supercavitating  model  moving  on  inertia; 

- the  basic  dimensions  of  unsteady  cavity  (length  and  mid  - section  diameter)  in  any  point  of  the 
distance. 

The  SCAV  program  (SuperCavity,  Version  1.3)  is  intended  to  automatize  designing  the  high-speed 
supercavitating  models. 

At  experimental  investigations  of  high-speed  motion  in  water,  the  model  shape  must  satisfy 
requirements  of  inscribing  into  a natural  cavity  and  structural  stability  of  the  supercavitating  motion 
regime.  If  we  suppose  that  the  model  moves  rectilinearly  under  action  the  cavitation  drag  only, 
then  the  process  is  described  by  the  dimensionless  equation  of  model  motion  on  inertia 


dV 

dt 


1 pwDncx0 


Fr 


8m 


2 — 2//  _ / \ 

V +Patm+ J-PcV) 

Fr2 
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where  m is  the  mass  of  model;  cxQ  is  the  cavitation  drag  coefficient  when  (7  = 0,  patm  is  the 
atmospheric  pressure;  H is  the  immersion  depth. 

Fig.  13,  a presents  a copy  of  PC-screen  after  calculation  of  the  model  location  in  the  cavity.  The 
model  shape  and  starting  data  correspond  to  the  real  experiment.  The  vertical  scale  of  the  image  is 
increased  for  clearness. 

In  Fig.  13,  b,  graphs  of  changing  the  cavity  length  Lc  and  model  velocity  V in  dependence  on  passed 
distance  x are  shown.  The  calculation  was  executed  for  titanic  model  at  Dn  = 1.5  mm,  m = 49.0  g, 
V0  = 1000  m/s.  It  is  supposed  that  the  model  enters  in  water  in  point  x = 0 and  then  moves  on  inertia. 
The  initial  rectilinear  part  of  the  graph  Lc  (x)  corresponds  to  the  interval  of  the  supercavity  formation. 
Quasistationary  magnitudes  of  Lc  (x)  are  plotted  by  dashed  lines  for  comparison. 


4.1.  Inertial  motion  with  propeller 


Motion  of  the  model  with  propeller  has  two  features.  First,  thrust  Tp  acts  on  the  model.  It 
compensates  the  drag.  Second,  the  part  of  the  model  material  with  density  ph  must  be  changed  by  fuel 

with  density  in  4 times  lower.  Therefore,  the  model  mass  with  propeller,  especially  in  the  path  end,  is 
lower  than  the  solid  model  mass. 

The  missile  thrust  can  be  determined  from  the  relation  behavior 

dm  , 

(30' 

where  T is  the  propeller  thrust;  m/  is  the  fuel  mass;  Q f is  the  specific  impulse  of  the  fuel.  Usually,  its 
value  is  expressed  in  kgf/(kg/sec)  or  lbt  /(lb/sec).  Numerical  magnitudes  are  equal  in  the  both  cases. 

Typical  magnitudes  of  Qf  are  in  range  from  200  to  500  kgf/(kg/sec).  Passing  onto  the  system  SI,  we 
obtain 


9.810, 


n-  sec 


9.81  Qf 


kg  ■ m sec 
sec2  kg 


m 


9-810, — 


sec 


Impulse  of  force  Tp  during  the  motion  time  is 


Kf  = Qfmf  . (31) 

If  the  propeller  chamber  volume  is  equal  to  Wc , then  the  model  mass  decrease  is 
A m0  =(pb—  p f )WC  in  beginning  of  the  path; 


A m0  = phWc  *n  the  path  end. 
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As  a result,  the  initial  impulse  is 


Kx  = [m0  - (pb  - p f )WC  ]V0 . (32) 

During  the  propeller  work  the  impulse  K , will  be  obtained  by  Eq.  (31),  and  general  available  impulse 
will  be 


K0  = Wc[pfQf-(ph-pf)V0]. 


The  expression  in  brackets  vanishes  at 


Pf 


Pt  ~ P, 


(33) 


When  V0  < Vcr , the  initial  impulse  increases,  and  we  can  hope  that  the  range  increases.  On  the  contrary, 
when  V0  > Vcr , it  decreases,  and  we  can  wait  for  the  model  with  propeller  passes  smaller  distance  than 
solid  one. 


4.2  Influence  of  a missile  propeller  on  the  model  motion  range 

Models  hi  - h4  with  length  250  mm  and  with  length  150  mm  were  chosen  as  base-models  to 
research  influence  of  the  missile  propeller  on  the  model  motion  range  (Tablel). 

Table  1 

Characteristics  of  models  of  a series  h 


N 

L 

Rn 

Re 

m„ 

me 

Wf 

Mf 

A m 

cm 

cm 

cm 

g 

g 

3 

cm 

G 

% 

% 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

hi 

25 

0.05 

0.7 

156 

115 

5.3 

98 

27 

6.3 

h2 

25 

0.15 

1.525 

749 

578 

21.0 

39 

22 

5.2 

h3 

15 

0.15 

1.225 

277 

227 

6.5 

12 

18 

4.3 

h4 

15 

0.15 

1.225 

277 

218 

7.5 

14 

21 

5.0 

Note:  m0  is  the  solid  model  mass;  me  is  the  mass  of  model  with  chamber;  W f is  the  chamber  volume; 
m f is  the  fuel  mass. 

It  was  accepted  that  the  chambers  having  truncated  cone  shape  were  made  in  models  to  dispose  the 
fuel.  A diameter  of  the  frontal  base  of  the  chamber  was  accepted  equal  to  the  half-diameter  of  the 
model  in  the  place  of  the  chamber  base  disposition.  A diameter  of  the  back  base  was  accepted  equal  to 
the  half-diameter  of  the  tail  part  of  the  model. 

As  a result,  four  models  were  developed.  They  are  designated  as  models  of  a series  h.  These  models 
characteristics  are  given  in  Tablel.  Their  schemes  are  shown  in  Fig.  14.  Models  h3  and  h4  differs  by  the 
chamber  dimension  only. 
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The  calculations  were  carried  out  by  using  the  program  CDILM.  We  varied  the  starting  velocity  of  the 
model  (1200,  1000,  800  and  600  m/s)  and  specific  impulse  of  the  fuel  (200,  400  and  600  kgf  /(kg/sec)). 
The  fuel  density  was  accepted  equal  to  1.85  g/cm3  in  all  the  cases. 

Calculations  for  the  case,  when  the  fuel  density  was  equal  to  the  model  material  density  (7.85  g/cm3 ), 
and  its  specific  thrust  was  equal  to  zero,  were  performed  for  comparison.  This  case  corresponds  to  the 
case  of  the  solid  model  without  propeller. 

It  is  accepted  in  the  calculations  that  the  combustion  occurs  at  constant  velocity.  All  the  fuel  had 
burnt  during  the  motion  and  the  fuel  combustion  continued  1/2  sec  to  1/12  sec. 

Thus,  the  model  was  passing  the  greatest  part  of  the  path  at  the  operating  propeller. 

One  can  see  changing  the  model  material  on  the  fuel  does  not  give  an  gain  in  the  range  for  any  model 
at  specific  impulse  0,  = 200.  The  most  effective  fuel  gives  some  gain  in  the  range  from  30%  at  high 
starting  velocity  to  150  % at  low  one  (Fig.  15). 


5.  Interaction  between  the  propeller  and  supercavity 


Problems  on  interaction  between  the  propeller  and  a supercavity  are  weakly  studied  direction  of  the 
hydrodynamics  and  enough  complex  for  theoretical  and  experimental  investigation. 

The  similar  problem  on  interaction  between  the  propeller  and  ship  is  known  in  shipbuilding  [16].  This 
interaction  may  be  both  positive: 

- the  system  efficiency  increasing  at  creation  of  the  additional  thrust  on  the  hull  by  means  of  the  pressure 
redistribution  and  negative; 

- the  entraining  force  increase  on  the  hull  and  the  interaction  efficiency  reduction  of  the  system  “hull- 
propeller”  in  dependence  on  construction  of  the  propeller  and  hull,  their  mutual  disposition  and  the 
propeller  work  regime. 

In  this  case  the  propeller  can  influence  on  the  cavity  shape  directly  by  changing  its  closure  conditions, 
free  surface  structure  and  gas-leakage  from  the  supercavity  at  its  presence  on  the  hull. 

The  following  types  should  be  noted  among  the  constructive  features  of  the  propellers: 

1.  Gas  jet. 

1.1.  Jet  with  condensable  phase. 

2.  Water  jet. 

3.  Screw. 

3.1.  Supercavitating  screw. 

4.  Turboprop. 


The  following  configurations  are  noted  among  the  configurations  of  propellers  about  the  supercavity: 
1 . The  frontal  disposition  of  the  propeller  or  the 
towing  propeller  located  in  the  nose  in  front  of  the  supercavity. 

2.  Tail  disposition  of  the  propeller  or  the  pushing  propeller  past  the  supercavity. 
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5.1.  Description  of  parameters  governing  the  dynamics  of  cavity  entrainment 


The  supercavity  as  the  flow  form  created  by  free  streamlines  is  very  sensitive  to  the  external  pressure 
fields,  including  the  pressure  fields  formed  by  the  propeller. 

As  is  known  from  the  hydromechanics  [16],  the  propellers  always  forms  jets  of  gas  or  fluid  modeled 
by  distributed  sources  and  dipoles. 

The  scheme  with  ideal  propeller  [16]  (Fig.  16)  gives  a typical  picture  of  change  of  the  flow 
parameters.  In  this  case  the  additional  fluid  is  not  injected  compared  to  the  case  of  the  jet  (Fig.  16),  and 
the  working  section  of  the  propeller  may  e modeled  by  a surface  of  dipoles  oriented  so  that  they  pumps 
the  fluid  from  the  zone  in  front  of  the  propeller  to  the  zone  past  the  propeller.  Thus,  the  pressure 
difference  from  — A P to  + A P is  created  in  the  working  section. 

It  is  clear  that  in  the  case  of  interaction  with  the  supercavity  the  frontal  zone  with  reduced  pressure 
— APwill  entrain  the  cavity,  and  the  cavity  will  be  displaced  from  the  zone  with  increased  pressure 
+ A P . It  is  necessary  to  consider  their  joint  influence  in  the  real  case  when  the  flow  is  formed  both  by 
the  hull  and  the  propeller. 

Since  the  propeller  must  only  compensate  the  impulse  loss  caused  by  the  hydrodynamic  drag  of  the 
hull  in  the  free  stream,  then  main  dimensionless  criterion  influencing  on  behavior  of  the  flow  around  the 
propeller  with  taking  the  hull  into  account 


T=^~,  (34) 

K 

where  T is  the  propeller  thrust,  R is  the  hull  drag.  In  this  case  it  is  necessary  to  distinguish  three 
regimes. 

When  T > 1 , it  is  the  acceleration  regime.  The  thrust  is  higher  than  the  hydrodynamic  drag  (the 
wake  past  the  body  has  the  positive  impulse); 

When  T — 1 , it  is  the  steady  motion  regime  (the  wake  without  impulse); 

When  T < 1 , it  is  the  deceleration  regime)  the  wake  past  the  body  has  the  negative  impulse). 

This  criterion  T is  the  most  considerable,  because  the  main  consumption  of  the  impulse  and  energy 
is  included  in  the  propeller  thrust.  This  results  in  the  basic  change  in  the  environmental  fluid. 

It  is  necessary  together  with  this  (34)  to  take  into  account  the  local  flow  parameters  caused  by 
propeller  directly  in  the  cavity  closure  zone. 


5.2.  Cavity  closure  in  the  gas  jet  zone  of  the  propeller 


• The  gas  jet  is  source  of  gas-supply  into  the  supercavity 


At  certain  conditions  of  gas  jet  outflow,  when  the  jet  outflow  velocity  F.  is  lower  than  the  critical 
velocity  VjCR 


- Vi  ~ 

V,  = — < V 
' V 


JCR  > 


(35) 


the  gas  part  may  return  from  the  jet  boundary  layer  into  the  cavity  and  be  the  additional  gas-supply. 


17-16 


• The  gas  jet  is  the  ejector  (sink)  sucking  out  the  gas  from  the  cavity 

This  regime  may  be  observed  when  Vj  exceeds  the  critical  velocity  VjCR 


- Vi  ~ 

Vj=J->VjCR, 


(36) 


The  criterion  F depends  on  the  place  and  conditions  of  the  cavity  closure.  So  the  criterion  V j (35), 
(36)  is  the  second  parameter  governing  the  dynamics  of  supercavity. 


5.3.  Interaction  between  the  cavity  and  propeller.  Tail  location  of  the  propeller 

Experiments  on  study  of  the  propeller  influence  on  the  cavity  were  carried  out  at  the  IHM  on  the 
special  installation  (Figs.  17,  18).  On  this  installation  the  supercavity  was  formed  by  means  of  the  ring 
cavitator  on  the  hull  with  diameter  Dm  = 52  mm  which  is  in  front  of  the  propeller  with  blades. 

Pictures  of  flow  around  the  model  at  different  regimes  of  the  propeller  work  at  continuous  gas-supply 
are  shown  in  photos  (Figs.  19,  20,  21). 


5.4.  Frontal  location  of  the  screw  propeller 


Investigation  of  the  propeller  screw  influence  on  the  cavity  was  carried  out  on  the  special  installation 
permitting  to  create  the  artificial  cavity  around  the  propeller  tail  cone  at  different  regimes  of  the  propeller 
screw  work.  A ratio  of  the  screw  and  cavity  diameter  was  Dp  = (2  3 )Dm  . The  propeller  screw  is 

located  in  the  nose  part  in  front  of  the  cavity.  It  is  established  in  virtue  of  the  fulfilled  researches  that  the 
working  propeller  screw  considerably  changes  the  cavitation  characteristics  of  the  tail  cone  [ 15]. 


Results  of  the  cavity  testing  for  the  working  propeller  screw  are  given  in  the  form  of  dependencies  of 


the  relative  parameters  of  the  artificial  cavity  as  function  of  the  screw  load  coefficient  K 


8 Kt 
7tX{ 


, where 


Kx  is  the  thrust  coefficient,  and  A;.  is  the  relative  screw  tread. 


Changes  of  the  cavity  length  L=Lp/L0,  cavitation  number  <7  = o / <7„  and  air  rate  coefficient 
Q =Qp/Q0,  where  L0 , a and  Q0  are  the  isolated  cavity  parameters  at  absence  of  the  propeller  screws 

but  at  other  the  same  conditions,  are  given  in  Figs.  22,  23,  24,  respectively.  As  is  shown  in  the  graphs, 
the  cavity  length  considerably  decreases  when  the  propeller  screw  is  located  in  front  of  the  cavity.  It  is 
about  30%  of  the  isolated  cavity  length  even  at  zero  load  of  the  screw  and  the  same  air  rate.  In  this  case 
the  cavitation  number  increases  but  considerably  slower  than  the  cavity  length  changes.  This  is 
explained  to  that  the  propeller  screw  not  only  increases  the  pressure  downstream  and  strongly  swirls  the 
stream.  As  a result,  the  cavity  boundary  is  destroyed.  To  reach  the  same  cavity  length  at  the  working 
propeller  screw  it  is  necessary  to  increase  in  some  times  the  air  rate.  So,  when  K = 0 , the  consumable 

air  rate  must  exceed  in  2-3  times  the  corresponding  value  of  Q0  for  the  isolated  cavity,  and  when 
K = 0.3 , it  is  necessary  to  increase  in  5 times  the  air  rate  (Fig.  24). 
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Properties  of  different  types  of  the  fuel 


Fuel 

Chemical 

Formula 

Density, 

g/cm3 

Molecular 

weight 

Temperature  of 
melting 

Temperature  of 
boiling 

Reaction 

Molecular 
weight  of  reaction 
product 

Quantity  of  evolved 
hydrogen , 

1/g. 

Specific  heat  of 
reaction, 

Kcal/1 

Theoretical  velocity 
of  flowing  out,  m/s 

Theoretical 
specific  impulse 

Aluminum 

ES 

2.7 

26.98 

658 

2500 

2A1+  1.502=  A1203 

101.94 

— 

3730 

5590 

570 

Aluminum 

E9 

26.98 

658 

2500 

A1  + 3H43  = Al(OHh  + 1 .5H, 

81.03 

1.25 

4150 

5900 

601 

Magnesium 

Mg 

1.74 

24.32 

657 

1102 

Mg  + 0.5O2  = MgO 

40.32 

- 

3623 

5510 

562 

ES£9 

1.74 

24.32 

657 

1102 

Mg  + 2H,0  = Mg(OH),  + H, 

60.32 

0.924 

3800 

5640 

575 

Lithium 

Em 

0.534 

6.94 

179 

1372 

2 Li  + 2H,0  = 2LiOH  + H2 

49.88 

1.62 

5340 

6700 

683 

Fig.  1 1 
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